Accepted by ApJ on Nov 11, 2009 



Galactic Cepheids with Spitzer: I. Leavitt Law and Colors 
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ABSTRACT 

Classical Cepheid variable stars have been important indicators of extragalac- 
tic distance and Galactic evolution for over a century. The Spitzer Space Tele- 
scope has opened the possibility of extending the study of Cepheids into the mid- 
and far-infrared, where interstellar extinction is reduced. We have obtained pho- 
tometry from images of a sample of Galactic Cepheids with the IRAC and MIPS 
instruments on Spitzer. Here we present the first mid-infrared period-luminosity 
relations for Classical Cepheids in the Galaxy, and the first ever Cepheid period- 
luminosity relations at 24 and 70 /im. We compare these relations with theoretical 
predictions, and with period-luminosity relations obtained in recent studies of 
the Large Magellanic Cloud. We find a significant period-color relation for the 
[3.6] — [8.0] IRAC color. Other mid-infrared colors for both Cepheids and non- 
variable supergiants are strongly affected by variable molecular spectral features, 
in particular deep CO absorption bands. We do not find strong evidence for mid- 
infrared excess caused by warm (~ 500 K) circumstellar dust. We discuss the 
possibility that recent detections with near-infrared interferometers of circumstel- 
lar shells around 5 Cep, £ Car, Polaris, Y Oph and RS Pup may be a signature 
of shocked gas emission in a dust-poor wind associated to pulsation-driven mass 
loss. 

Subject headings: Cepheids — distance scale — infrared: stars — stars: mass 
loss 
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Introduction 



Although it was over 100 years ago tha t Henrietta Leavitt discovered the Cepheid 
Period-Luminosity relation (PL. iLeavittl Il908l ). or "Leavitt Law", few tools in astronomy 
have such enduring importance. Classical Cepheid variable stars are fundamental calibrators 
of the extragalactic distance scale and in addition their observed properties are a benchmark 
for stellar evolution models of intermediate mass stars. It is now widely accepted that 
distances derived from SN la are still significantly influenced by the classical Cepheid cal- 
ibration: iRiess et al.l (120051 ) found a change of 15% in the value of Hq when only modern, 
high-quality SN la data and HST Cepheids were used. Recent works have extended optical 
and near-infrared PL and period-color (PC) relations to the mid-infrared (mid-IR), where 
there is le ss interstellar extinction, with observations o f Large Mage l lanic Cloud (LMC) 



Cephe ids (jFreedman et al.l 120081 ; iNgeow fc Kanbur 



2009al ) with the Spitzer Space Telescope ( jWerner et al.l 120041 ) 



2008 



Ngeow et al.l 1200a Madore et al. 



Evolutionary and pulsation properties of intermediate-mass stars in the core He-burning 
phase, like Cepheids, play a crucial role in several long-standing astrophysical problems. 
They are transition objects between stellar structures ending up their evolution either as 
white dwarfs or as core collapse supernovae. Therefore, they are not only the most popular 
primary distance indicators, but are also crucial to understanding the chemical evolution of 
stellar systems hosting a substantial fraction of young stars, e.g the Galactic disk and the 
dwarf irregular galaxies in the Local Group and in the Local Volume (d < 10 Mpc). 



Although these obj ects are fundamental for stellar evolutio n (jBono et al.ll2000l ; iBeaulieu et al. 



20011 ) . stella r pulsation (|Bono et al 



1999; Ma rconi et al.ll2005h . and Galactic chemical evolu- 



tion models ( Pedicelli et al. 2009 ; Spitoni et al. 20091 ). current predictions are still hampered 



by several problems. The most outstanding issue is the "Cepheid mass discrepancy" between 
pulsation masses of classical Ceph eids and thei r evolu tionary masses. Evidence was brought 
forward more than 30 years ago by iFricke et al.l (119721 ) who found that pulsation masses were 
from 1. 5 to 2 times smaller than the evolutionary masses. This conundrum was partially 



solved (IMoskali k et a 



Project (jSeaton et al 



1992J) by the new set s of radiative opacities r eleased by the Opacity 



19941 ) and by OPAL (IRogers fc Iglesia s 1992) . However, severa l re- 



cent investigations focused o n Galactic Cepheids (Bono et al 



Magellanic Cloud Cepheids (IBeaulieu et al.l 12001 



2001 



Bono et al 



Caputo et al.ll2005j) and 



2002 



Keller fc Wood! \200& \ 



suggest that such a di screpancy still am ounts to 10-15%. Measured masses of Galactic 
binary Cepheids (e.g. lEvans et al.l 120081 ) are also smaller than predicted by evolutionary 
models neglecting core convective overshooting during central hydrogen burning phases. 

The relative importance of the main factors affecting Cepheid model mass estimates 
(extra-mixing, rotation, radiative opacity, mass loss and binarity) is still debated. Even 
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t houg h commonly used semi-empirical relations (IReimers et al.lll975c iDeJager &: Nieuwenhuijzen 
19971 ) do not predict enough mass loss to solve the Cepheid mass discrepancy problem, mass 
loss may indeed be the key culprit among the physical mechanisms suggested to explain the 



mass discrepancy problem. The semi-empirical mass loss relation derived by lReimers et al. 



( jl975l ) is clearly inade quate to correc tly estimate the mass loss for certain evolutionary phases 
of giant stars (see e.g. IWillsonll2000l ). A plausible increase in the typically adopted Reimers 
wind free parameter (0.2 < r\ < 0.4) does not account for t he entire range in mass covered by 



clust e r horizontal branch (cen tral helium burning) stars (lYong et al.ll2000l ; ICastellani et al. 



20051 ; ISerenelli &: Weissl 120051 ) . There is no good reason why a similar discrepancy should 
not affect intermediate-mass stars. 

Mass loss is often betrayed by the presence of dust shells, detectable as infrared excesses, 
or by stellar winds which produce blue-shifted absorption dips in the ultraviolet. Empirical 
estimates of mass loss rates based on infrared (IRAS) and ultraviolet (IUE spectra) obser- 
vations for a large sample of G alactic Cepheids suggest mass- loss rates ranging from 10~ 10 
to 10~ 7 M Q yr _1 (IDeasyl 119881 ) . However, evidence for mass loss rates high enough to af- 
fect evolution is very rar e, and it is not clear tha t mass loss is a wide-s pread phenomenon. 
McAlary fc Welch! (119861 ) used IRAS photometry (jBeichman et al.lll985l ) and found evidence 
of very cool dust {T d < 50 K) around two classical Cepheids (RS Pup and SU Cas) known 
to be associated with reflection nebulae. Due to the very large IRAS beamsize (as large as 
~ 5 arcmin), however, it was very difficult to separate local dust emission from "Galactic 
cirrus" background emiss i on. More recently, K band nea r-IR interferometric observations 
jMerand et alibood . 120071 : iKervella et alihood . liooi . bood ) have detected circumstellar emis- 
sion around five Classical Cepheids: 5 Cep, £ Car, Polaris, Y Oph and RS Pup. While the 
nature of the material responsible for this emission remains mysterious, this is a tantalizing 
suggestion that mass loss activit y may be pres e nt ar ound these nearby Cepheids. A similar 
conclusion was also reache d by iNeilson et al.l (120091) . that using OGLE (Optical Gravita- 
tional Lensing Experiment; Udalski et al.1l999 ) and Spitzer data for Magellanic Cepheids, 
found evidence of a wide range of mass loss rates. 

To investigate this possibility, we have obtained Spitzer observations of a sample of 
Galactic Classi cal Cepheids. All stars were observed with both the Spitzer Infrared Array 
Camera (IRAC.lFazio et al.ll20041) and the Multiband Infrared Photometer for Spitzer (MIPS, 
Rieke et all [20041 ). The aim of the IRAC observations was to characterize Galactic Cepheid 
colors in the mid-IR, where the emission from the photosphere is still dominant, and search 
for infrared excess related to warm (~ 500 K) circumstellar dust. The MIPS observations 
were intended to investigate the presence of extended emission from cool (< 100 K) dust, 
taking advantage of the higher angular resolution of Spitzer (5 arcsec at 24 /xm). In this 
paper we present the photometry of our Cepheid sample in the IRAC and MIPS bands, 
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discussing their PL and period-color (PC) relations. The results of our search for extended 
emission in IRAC and MIPS images will be presented in a separate paper (Barmby et al. in 
preparation). 

The criteria for our sample selection are laid out in section [2J, and the observations 
are described in section [31 The techniques adopted to measure the source photometry in 
all IRAC and MIPS bands are discussed in detail in section [U In section [5] we derive the 
Leavitt Law and PC relations in the IRAC and MIPS bands, and we compare them with 
similar relations obtained for the LMC. We study the intrinsic mid-IR colors of Cepheids 
in section [61 where we set limits on infrared excess at Spitzer wavelengths. Our results are 
discussed in section [7] and summarized in section [HI 



2. Sample Selection 



Our sample stars are listed in Table [H They include 29 bright (K < 5.1) nearby 
Cepheids, as well as 3 non-variable supergiants and one red giant added for comparison. 
The sample was selected to cover a range of several characteristics which may influence mass 
loss, or enable us to understand it better. Several of the Cepheids are in clusters, three 
are first-overtone pulsators, and eight are in multiple star systems. The range of periods is 
about a factor of ten, and several Bump Cepheids (9 < P < 12 days) are included so that 
we could investigate whether these objects present a pec uliar mass loss rate. Our sample 



includes all 9 Cepheids wi th measured angular diameters (INordgren et al.ll2000l ; lLane et al. 



20021 ; iKervella et al.ll2004l ) at the time of the Spitzer proposal submission , as well as 5 stars 
(Polaris, S Mus, W Sgr, V350 Sgr, an d FF Aql) with observed masses ([Evans et al.l 119971 
20061 . 120081 . 120091 : iBenedict et al1l2007h . V350 Sgr, S Mus, RS Pup, and f3 Dor are the only 
Cepheids with 60 /zm excesses observed by IRAS, so they can be used to check the earlier 
results. Our sample also includes Y Oph, a long period Cepheid with a surprisingly small 
amplitude and placid velocity structure (based on line profiles). Its unusual characteristics 
make it a candidate for an unusual history (either mass loss, coalescence, peculiar chemical 
composition, evolutionary status). 

The stars were also selected for the availability of accurate distance determinations (with 
the exception of the cluster star V636 Sco whose parallax is very uncertain). For 7 Cepheids 
we ha ve adopted the HST Fine Guidance Sensor parallaxes me asured by 



Benedict et al. 



(120071 ) . while we have adopted the Hipparcos distance for Polaris (Ivan Leeuwen et al. I l2007h . 



For the other stars we have adopted distances derived with the infrared surface brightness 
(hereafter IRSB) technique (see e.g. iFouque &: Gierenlll9971 ). While the IRSB method can 
provide very precise parallaxes when the pulsational velocity and lightcurve of a pulsating 
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star is known, the conversion between the observable radial velocity v r into pulsational 
velocity v p depends on the so-called projection factor p = v p /v r . The p-factor, despite its 
name, does not depend exclusively on t he geometry of the s tar, but also on its physical 
properties, including the period (see e.g. iNardetto et al.l 120071 ). Several p(P) relations have 
been proposed, and agreement on whi ch one better describes the true dependence of the 
p-factor from period is missi ng (see e.g. iRomaniello et al.ll2008l). Until 2007 several catalogs 



of Cepheid IRSB distances ( jStorm et al. 



2004 



Groenewegen et al.ll2004l; [Barnes et al.l 12005 



among others) used the relation p = 1.39 — 0.03 log P (iGieren et al.lll993f) characterized by 
a weak dependence from the period. More recently INardetto et al. I J2007h p ropos ed a new 
"steeper" relation p = 1.366 — 0.075 logP, which is used in the IFouque et al.l (120071 ) catalog. 



In Table [T] we provide two separate columns for the di stance of our target s: the first 
column mainly lists the "new" distance determination from IFouque et al. (b007h . while the 
second lists the "old" IRS B distances. In both ca ses, whenever available, we adopt the 



Benedict et al.l ( 120071 ) and Ivan Leeuwen et al.l (120071 ) distances, which are independent from 
the p-factor. For three stars (DT Cyg, V350 Sgr and U Aql) we could not find their IRSB 
distance estimated using the "old" p(P) relation. For completeness we derive the "old" 
distance of these stars (phis V636 Sco) using the method described in Fouque k. Gieren 
( 119971 ). These distances, and the parameters used to derive them, are listed in Table [2j 
V and E( B — V) are derive d from the Dunlap O bservatory Databas e of G ala ctic Classical 



Cephe ids (jFernie et al.lll995l ). K magnitudes from Ivan Leeuwen et al.l (120071 1 or lWelch et al. 
jl984h . while R and AR from iMoskalik fc Grovnval (b()05h . 



All stars except 3 (DT Cyg, S Z Tau and Polaris) are listed as fundamental mode pul- 
sators, according to IFouque et al.l (120071). FF Aql is here assumed to pulsate in the funda- 
mental mode, as proposed bv lBenedict et al.l (120071 ) , desp ite having been previously classified 



as a 



;ai moae, as proposed Dy Benedict et al. (zuuvj) , aesp i 
first overtone Cepheid by Feast fc Catchpolel ( 19971 ). 



3. Observations 

The observations were executed between July 19, 2006 and October 28, 2007 as part of 
the Spitzer Space Telescope Cycle 3 General Observer program with PID 30666. Each star 
was observed in IRAC 3.6, 4.5, 5.8 and 8.0 /an and MIPS 24 and 70 /jm bands. The MIPS 
160 jum band was not used: the known near-infrared leak in this filter would have resulted 
in strong contamination by starlight. 

The IRAC observations were designed to provide, whenever possible, high S/N unsatu- 
rated images at all wavelengths well within the linearity regime of the detectors, using the 
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IRAC subarray mode with 0.02 or 0.4 sec frame times. For targets brighter than K ~ 2.5 
even the shortest subarray frame time was too long to prevent saturation. For this reason, 
6 stars were observed using IRAC full frame stellar mode (frame time of 0.4 sec at 3.6 and 
4.5 /im and 2 sec at 5.8 and 8.0 jum), with the intent to produce heavily saturated images 
suitable for PSF-fitting photometry on high S/N unsaturated "wings" and diffraction spikes. 
Four other targets were observed in both subarray and full frame stellar mode, to provide 
photometric cross-calibration between the two modes of observations and to check for linear- 
ity. Polaris was not observed with IRAC as part of this program, since archival data obtained 
during Cycle 1 (Spitzer program 19) was available. Each star was dithered on the array in 4 
(subarray) or 5 (full frame) positions, using a small scale Gaussian offsets pattern, in order 
to facilitate transient outlier and bad pixel removals, and produce a better spatial sampling 
of the images. The total time on source, for each band, was 25.6 sec for the subarray observa- 
tions. Full frame observations had 2.0 sec (3.6 and 4.5 /im) and 10 sec (5.8 and 8.0 /im) total 
integration times. The IRAC data were reduced starting from the Basic Calibrated Data 
(BCDs) generated by the Spitzer IRAC pipeline versions S14.4.0 through S16.1.0. Mosaic 
images with a pixel scale of 0.8627 arcsec /pix for each source in each band w ere then created 



using the IRACproc post-BCD software (ISchuster. Marengo fc Pattenll2006l ). 



MIPS observations at both 24 /im and 70 /im were obtained with the Photometry 
Astronomical Observation Template. The 24 /im observations used 1 cycle of 3 sec frames 
with the standard 16 point dither pattern. The MIPS 70 /im observations used the default 
pixel scale and small field size, with 1 to 12 cycles of 3 or 10 sec frames (depending on source 
and background brightness) with the standard dither pattern. The cluster Cepheids S Nor, 
U Sgr, and V Cen had slightly different 70 /im observations: 2 adjacent fields of view were 
observed, to increase the sky coverage to 5 x 5'. The total on-source integration times were 
48.2 sec at 24 /im and from 37.7 to 1,300.2 sec at 70 /im. The MIPS data, generated by 
the Spitzer pipeline version S16.1, was retrieved from the Archive. No post-processing was 
performed on most of the data: post-BCD mosaics were used for the 24 /im data and most 
of the 70 (Ltm data (filtered versions were used). For the three cluster stars (S Nor, SZ Tau, 
U Sgr), the 70 /im data was remosaiced using MOPEX to combine the data for the two 
fields of view. The pipeline mosaics for two other stars (RS Pup, GH Lup) showed negative 
sidelobes; for these stars the BCDs were time- and column filtered and the mosaics re-made. 



4. Source Photometry 



We measured the IRAC Vega magnitude of all sources observed in subarray mode using 
aperture photometry. We derived photometric zero points from the IRAC absolute calibra- 
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tion factors (FLUXCONV) in the BCD headers (listed in Table [3]). For most sources we 
used an aperture of 12.2 arcsec radius (~ 10 pixels in the IRAC array pixel scale) with 
background level determined from a sky annulus with 12.2 and 24.4 arcsec inner and outer 
radii. This aperture is the standard used for the IRAC point source calibration, and as such 
corresponds to an aperture correction equal to one. Three stars (AQ Pup, VY Car and 
FF Aql, see Figured]) had one or more background sources falling within this aperture, and 
required using a smaller 6.1 arcsec radius aperture (~ 5 IRAC pixels). For these sources we 
derived our own aperture corrections, listed in Table [3J using the other stars as templates. 
This determination of the aperture correction is more reliable than using the values listed in 
the IRAC Data Handbook^ (2006), which were obtained from images with different spatial 
sampling and dither pattern. We have corrected our subarray point sou rce photom e try fo r 



the position-dependent geometric and gain distortion factors provided by lHora et al.l (120081 ). 
These factors (f CO rr) are listed in Table [3] for the subarray central position. The photometric 
uncertainty was estimated adding in quadrature the sky RMS variation and source pho- 
ton noise. At 3.6 and 4.5 jum, an extra 1% uncertainty due to uncharacterized pixel-phase 
variations was also added (see the IRAC Data Handbook, 2006, p. 50). 

The IRAC photometry of the stars observed in full frame stellar mode was obtained by 
fitting the unsaturated "wings" and diffraction spikes of the mosaic images with the IRAC 
extended Pixel Response Function (PRF). The extended PRF was constructed from multi- 
epoch deep images of Sirius, Vega (2 epochs), Fomalhaut (2 epochs), e Eridani (2 epochs) 
and e Indi, obtained in a separate Guaranteed Time Observation program (PID 90) and 
IRAC calibration observations. The individual observations were reduced with IRACproc, 
and the mosaic of each star was masked to remove any pixel closer than 50% or 80% to 
the saturation level to preserve flux linearity. The individual images were then rescaled to 
one of the images of Vega and combined with a sigma-clipping algorithm in order to remove 
background sources. The final product is a high S/N representation of the extended features 
of the IRAC Point Spread Function (PSF) projected on the IRAC pixel grid (e.g. a PRF), 
and is ideally suited to derive PSF fitting photometry of bright saturated sources. Given 
that the PRF is normalized as a calibrated image of Vega, the best-fit scaling factor between 
a star image and the PRF is equal to the flux ratio of the star with Vega, and can be 
directly converted into Vega magnitudes without relying on the IRAC absolute calibration. 
Magnitudes derived with this procedure are not affected by the position-dependent geometric 
and gain distortion factors, since the sources and the PSF stars have all been observed at the 
center of th e IRAC arrays. An ima ge of the IRAC PRF obtained with a similar procedure 



is shown in iMarengo et al.l (j2009al ). and the PRF files are available at the Spitzer Science 



J http: / /ssc. spitzer. caltech.edu/irac/dh/ 
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Center web sit^E 

The photometric error was estimated by bracketing the best subtraction with an over- 
and under-subtracted fit in which the PSF-subtraction residuals were of the same magnitude 
as the background noise. The fit was done using only the overlapping region where both 
the objects in our programs and the PRF were not saturated. This region is relatively 
narrow (width of ~ 20 - 30 pixels) at 3.6 and 4.5 jttm, due to the short total exposure time 
(2 sec) in these two bands. As a consequence, The PSF-fitting photometric errors tend to 
be larger at 3.6 and 4.5 fim than at 5.8 and 8.0 /zm. The dominant source of uncertainty 
for the PSF-fitting photometry is the "sampling noise" (grid pattern due to the coarse 
sampling of the IRAC PSF on the pixel grid) in the subtraction residuals. This coherent 
noise does not follow Poisson statistics. For this reason the photometric errors estimated for 
the PSF-fitting photometry should not be taken as la statistical errors (as in the case of the 
aperture photometry uncertainties), but rather as upper and lower limits in the photometric 
measurements. 

To test the agreement between aperture and PSF-fitting photometry we used the 4 stars 
observed in both modes. In general, the full frame and subarray mode observations of these 
stars have been executed at different epochs. As a c o nseque nce, based on Large Magellanic 



Cloud (LMC) Cepheid photometry by lMadore et al.l (l2009bl ). we should expect a variability 



in the IRAC bands as large as ±0.4 mag. Still, we can use these stars observed in both modes 
to assess the cross calibration of the aperture and PSF-fitting methods, although with this 
limitation. Table H] and Figure [2] show that we do not measure in our data any statistically 
significant shift between the subarray and the full frame photometry: the average difference 
between subarray and full frame magnitudes is much smaller than its RMS variation. The 
RMS variation (0.02, 0.07, 0.05 and 0.02 mags at 3.6, 4.5, 5.8 and 8.0 jttm respectively) 
is smaller than the average uncertainty in the PSF-fitting photometry estimated with the 
fit-bracketing method. For V Cen (the only case in which both subarray and full frame 
observations were executed at the same epoch), the difference between the aperture and 
PSF fitting photometry is 0.01, -0.04, -0.02 and -0.02 mags at 3.6, 4.5, 5.8 and 8.0 fjm 
respectively. This is within the uncertainty of the PSF fitting photometry of this star (0.27, 
0.16, 0.09 and 0.06 mags, respectively). At least in the case of V Cen, this test can be 
viewed as an independent confirmation of the IRAC absolute photometric calibration, since 
aperture photometry depends on the FLUXCONV parameters, while PSF-fitting with our 
PRF does not. 

Table [5] shows the final magnitudes of all sources. When available, we list the IRAC 



2 http: / /ssc. spitzer.caltech.edu/irac/psf.html 
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aperture photometry of unsaturated images (which has smaller nominal uncertainty) rather 
than the PSF-fitting photometry. The PSF-fitting photometry is instead listed for the sources 
observed only in IRAC full frame mode. 

Measurements of MIPS 24 /im flux densities followed standard procedures. Flux densi- 
ties for most objects were measured in a 20" radius apertures with background subtraction 
from sky annuli of 40"-50" and the recommended aperture correction of 1.13 applied. There 
are a few stars with apparent extended emission nearby; for these, flux densities were mea- 
sured in 7" radius apertures with the same sky annulus, and an aperture correction of 1.61. 
The flux densities measured for Polaris and £ Car (~ 3.7 Jy) are formally above the satu- 
ration limit for 3 second exposures, but both objects are below the "hard" saturation limit 
(4.1 Jy): the pipeline correctly replaced the saturated pixels with values from the 0.5 second 
exposures taken as part of the photometry AOR. Uncertainties for the 24 //m flux densities 
were computed using the standard IRAF/PHOT formula as: 



a 2 = (F/G + Aa 2 + A 2 a 2 /A 



l/Nl/2 



(1) 



where F is the flux in image units (MJy sr _1 ), G is the effective gain (the product of 
the detector gain, individual frame exposure time, and number of frames, divided by the 
FLUXCONV factor which converts from DN s -1 to MJy sr -1 ), A is the aperture area, 
A' is the background annulus area, and a is the standard deviation of image counts in 
the background area. T he values listed in Tabl e \5\ do not include the 2% uncertainty in 
the absolute calibration (jEngelbracht et al.l 120071 ). the dominant systematic error for these 
bright stars. 

Measurements of MIPS 70 /im flux densities also followed standard procedures. Flux 
densities were measured in 16" radius apertures fixed on the 24 /im position, with sky annuli 
of 18"-39", and an aperture correction of 2.04. There were no saturation issues in the 70 /im 
data. Uncertainties in the 70 fim flux densities were calculated using a similar equation 
to that for MIPS- 24, except that for 70 /im data Poisson noise is negligible. Following 
Carpenter et al.l (120081 ). the uncertainties are increased by extra multiplicative factors of 
^corr = 2.5 (noise correlation between pixels due to re-sampling during mosaicing), and 
r^sky = 1.5 (excess sky noise due to the data-tak ing procedure). The absolute calibration 
uncertainty at 70 /im is 5% ([Gordon et al.l 120071 ) . again not included in the data table. A 
number of stars were not detected at 70 /im: for these sources 3cx upper limits are recorded. 

Determining the effects of interstellar extinction on Spitzer IRAC and MIPS photometry 
is tricky, because of the difficulty of estimating the reddening in the Galaxy, c ompared to 



that for an external galaxy. Measurements of the Galactic extinction curve (ILutz et al. 
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19961 ; llndebetouw et al.ll2005t iFlaherty et al. 120071; iRoman-Zuniga et al.ll2007l) when scale d 



to IRAC wavelengths give A 5 ^ m /A B « 0.04 (IRieke fc Lebofskylll985l ; lFreedman et al.ll2008l ). 



We should thus expect reddening at mid-IR wavelengths for all our nearby Cepheids to be 
negligible, and for this reason we do not correct our data for reddening. 



5. Leavitt Law and Period-Color Relations 

The advantages of deriving a PL relation (the Leavitt Law) in the mid-IR (rather than in 
the optical or near-IR) are manifold: luminosity variation amplitudes, metallicity effects and 
interstellar extinction are all expected to be smaller at longer wavelengths. A well calibrated 
PL relation in the IRAC and MIPS bands is of great utility for the available Galactic and 
extragalactic stellar photometric catalogs obtained in the Spitzer legacy projects. 

The PC relation at optical wavelength provides important diagnostics for basic stellar 



parameters such as eff ective temperature and metallicity (see e.g. ISandage et al.ll2004l 12009 



Tammann et al. 20031 for a comparative study of the PC relations between the Galaxy, the 



LMC and the SMC). Extending the study of PC relations into the mid-IR offers the chance 
to study the effects of these parameters on different spectral features, at wavelengths where 
the interstellar reddening is less important. A firm understanding of the "intrinsic" color 
properties of classical Cepheids, and their dependence on stellar parameters and period, is 
also essential to search for the presence of "extrinsic" color excess that may be related to 
circumstellar dust emission. 



5.1. PL Relations 



An IRAC PL relation has already been obtained for Cepheids in the LMC (IFreedman et al 



20081 : iNgeow fc Kanbujjiool 



SAGE project ( jMeixner et al. 



2006h 



Ngeow et al.ll2009i ; iMadore et al.ll2009al ). using data from the 



and was used to derive the distance modulus for two 



galax ies in the Local Group, NGC 6822 (IMadore et al.ll2009bl ) and IC 1613 (IFreedman et al. 



20091 ) . Using our photometry in Table we can derive the same relation for a sample of 
Cepheids in our Galaxy. These are the same Cepheids that, by virtue of their independent 
distance determination, are at the base of the PL relation zero-point calibration. 

We derived the best-fit PL relations for all IRAC and MIPS data with a linear least 
square fit method, weighed on the uncertainty of each photometric point. For IRAC subarray 
data we use as weight the 3a photometric uncertainty, while for all other measurements we 
use the uncertainty quoted in Table [51 In particular, for IRAC full frame photometry we 



use the PSF-fit "bracketing" interval, which is similar to a 3a uncertainty. For MIPS 70 /im 
data we exclude data points for which only an upper flux limit is available. Table [6] lists the 
coefficients of the best-fit relations, in the form: 



M = a[\ogP -1.0]+ P (2) 

Note how the PL relations we have derived are based on single epoch measurements, and 
not on mean photometry that would have required multiple epoch observations. To include 
the first overtone Cepheids (DT Cyg, SZ Tau and Polaris) in the fit, we have converted their 
observed period Pi into a "fundamentalized" period P using the relation: 



Pi/Po = 0.716 - 0.027 log Pi 



(3) 



from iFeast fc Catchpold (119971 ) , derived in turn from lAlcock et al.l (119951 ) . We provide the 
best-fit PL relations for three different cases. First, using the "new" IRSB distance de- 
termination rel ying on the "s t eep" p(P) relation. Then, u s ing o nly the stars for which 
the astrometric Benedict et al.l (120071 ) and Ivan Leeuwen et al.l (120071 ) distances are available. 
Last, using the "old" distances based on the p(P) relation with weak period dependence. 
Note than in both "old" and "new" cases we still use the astrometric distances when avail- 
able. The best-fit PL relations for the "new" (solid lines) and astrometric (dashed lines) 
distances are shown in Figure [3J 



Table [6] also includes the K band relations from iFouque et al.l (120071 ) , Benedict et al. 



( 120071 ) and IStorm et al.l (12004 ) for the "new," astrometric, and "old" distances respectively. 
The table shows that the slopes and zero points of the three cases are consistent within their 
uncertainties. The astrometric PL relations have larger uncertainties because most of the 
stars in that sample have the less precise full frame PSF-fitting photometry, and because of 
the smaller size of the sample. We can however note some trends: the astrometric Leavitt 
laws are consistently less steep than the relations using IRSB distances. The "old" distances 
produce steeper relations than the "new" distances. A similar trend is present in the K band 
relations. The zero points are very similar for all cases, well within their uncertainties. 

Figure [3] also shows that FF Aql (data-point at logP = .73) closely matches the 



general trend. If instead, as suggested by lFeast fc Catchpold ( 119971 ) . the source was pulsating 
in the first overtone, its "fundamentalized" PL ratio would be a significant outlier. This 
result suppor ts the view that FF A ql is indeed behaving as a fundamental mode pulsator as 
suggested by Benedict et al.l (120071 ) . 



The dispersion around the best-fit PL relations is ~ 0.2 mag in the four IRAC bands 
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and MIPS 24 /zm, and ~ 0.4 mag at 70 /zm. This is comparable with the expected dispersion 
at Spitzer wavelengths, once our photometric and distance uncertainty is taken into account. 
The larger uncertainty in the best-fit parameters at 70 /zm is due to the larger photometric 
errors and the smaller number of sources for which a measured flux is available at this 
wavelength. 

The dispersion around the best-fit PL relations has correlated and uncorrelated com- 
ponents between bands. The correlated scatter is a consequence of Cepheid variability. The 
uncorrelated scatter <ta is instead due in part to the absolute magnitude uncertainties (ge- 
ometric sum of the photometric error and the uncertainty in the distance) and in part to 
intrinsic variations related to the "width" of the Cepheid instability strip (due to differences 
in metallicity and other stellar parameters). Figure H] (left column) shows that the correlated 
component of the fit residual has an amplitude of about —0.3 to +0.5 mag. We may assume 
this value to approximate the amplitude of the Cepheid light curves in the IRAC bands. 
While the weighted average residual is zero (by virtue of the least-square fit used to derive 
the PL relation), it is interesting to note that the IRAC residuals are asymmetrically dis- 
tributed. This may be a reflection of a prevalence of asymmetric light curves in our sample. 
While the amplitude of the correlat ed component of the P L dispersion in our data is similar 



to the one measured in the LMC by lFreedman et al.l (120081 ). they do not note any asymmetry 



in their sample (three times larger). Comparison between our measured a a (~ 0.04-0.05) 
and the absolute magnitude uncertainty (on average ~ 0.15 mag) shows that we do not have 
enough accuracy to measure the intrinsic scatter in the PL relation, and all the observed 
scatter can be explained by the uncertainty of our data. 

Our sample does not show any correlation between variability an d color residua l s in the 



IRAC bands (Figure IU right column), similar to the result obtained by lMadore et al.l (j2009al ) 
in the LMC. Again, we cannot test for color correlation between the IRAC and MIPS bands 
due to the different epochs at which the two sets of data were acquired. 



5.2. PC Relations 

Figure [5] shows the period vs. Spitzer color trends in our data. These diagrams are better 
suited to analyze color variations between the sample Cepheids than the color residuals 
in Figure HJ because they do not suffer from the larger uncertainties due to the distance 
estimates necessary to derive the PL fits. Even in the case of colors using only IRAC bands 
(taken simultaneously, and hence with each star at a given phase), the data show a scatter 
that is well above the photometric error. The only exception is the [3.6] — [8.0] color, for 
which it is possible to derive a weighted, least-square best-fit PC relation: 
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[3.6] - [8.0] = 0.039 [log(P) - 1.0] [±0.008] - 0.019[±0.011] (4) 

The only apparent outliers are some nearby Cepheids observed in full frame mode that 
have large photometric errors, and in particular Polaris which is consistently red in all IRAC 
colors by ~ 0.1 mag. This detection is above the PSF-fitting photometry uncertainty for 
this source, but should be taken with care, given that this observation was not planned to 
maximize the accuracy of the PSF fit, and may suffer from uncharacterized systematic errors. 
The large scatter in the [3.6] — [4.5] and [3.6] — [5.8] color (as large as ~ 0.1 mag) prevents 
the determination of a meaningful PC relation in these colors. This scatter (which tends to 
be larger for long period Cepheids) may be related to either some intrinsic property of the 
stellar emission (the presence of broad spectral features in their mid-IR spectrum) or some 
external cause (emission/absorption features of interstellar or circumstellar origin). These 
two possibilities will be examined in detail in the next section. 

Colors combining IRAC and MIPS bands have an additional scatter due to the different 
epochs at which the data were acquired. We can estimate this scatter to be at least ±0.4 mag, 
as determined in the previous section for the IRAC bands (assuming that the variability at 
MIPS wavelengths is negligible, which is not necessarily the case). To be detected, any 
excess in the [3.6] — [24] and [3.6] — [70] colors (or other combinations mixing IRAC and 
MIPS bands) must be larger than 0.4 mag. Figure [5] shows that none of our sources have 
excess in the [3.6] — [24] color above this threshold (dashed line). None of the sources with a 
70 jum measured flux has a [3.6] — [70] excess above 3a of the minimum expected variability 
scatter (dotted line), with only sources undetected at 70 /im above that line. Similar results 
are obtained by combining the MIPS 24 and 70 /j,m bands with the other IRAC channels. 

The [24] — [70] color does not suffer from multi-epoch scatter, but the lack of good 
photometric measurements for the 70 /im band prevent us from determining reliable colors 
for most of the sources. Only one target (SZ Tau) has a non-zero [24] — [70] color at more 
than 3a significance, which may indicate the presence of a true color excess from this star. 
The significance of this detection will be described in the next paper of this series (Barmby 
et al., in preparation). 

6. Cepheid Mid-IR Colors and Search for Infrared Excess 

Dusty mass loss in the wind of a Cepheid star would reveal itself in the form of an 
infrared excess caused by thermal radiation from the dust grains. IRAC and MIPS are 
particularly sensitive to this excess, as shown in the case of Asymptotic Giant Branch (AGB) 
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stars (IMarengo et al.ll2008l ). To isolate any excess from our data, we need however to first 
characterize the intrinsic colors of Cepheids at mid-IR bands. While a systematic study of 
the infrared spectral properties of Classical Cepheids is beyond the scope of this work, we can 
us e available mod e ls of t he 10 d period Cepheid £ Gem (one of our target stars) computed 



by lMarengo et al.l (120021 ) as a test case. 



Marengo et al.l (120021 ) developed a method to derive time- dependent hydrodynamic 
models of classical Cepheids by solving Local Thermo dynamical Equilibrium (LTE) plane 
parallel radiative transfer for a dynamic stellar atmosphere. The dynamic atmosphere, com- 
puted in spherical geometry and in non-LTE conditions, was in turn driven by a "pulsation 
piston" reproducing the correct radial profile of the star. While this procedure was not simul- 
taneously solving the hydrodynamic and radiative transfer solution for the star, it was still 
able to produce realistic time-dependent intensity spectra, then used to evalua te the temporal 



varia tions of the stellar limb darkening at optical and infrared wavelengths (IMarengo et al. 



20031 ). The method was applied to the case of ( Gem to estimate the effect of limb darkening 



variations on the star 's distance determination with the geometric Baade-Wesselink method 



(IMarengo et all 12004 ) . 



Figure [6] shows the hydrodynamic model spectra of ( Gem at 5 significant phases: 
4>l = (maximum visible luminosity), <pL = 0.32 (maximum radius), (pi = 0.49 (minimum 
visible luminosity), <pL = 0.63 (a phase at which a shock-wave is crossing the photosphere), 
<pL = 0.74 (minimum radius) and <pL = 0.85 (the phase at which the star was observed with 
IRAC). The thick line shows the hydrodynamic model spectra, while the thin line shows an 
equivalent hydrostatic equilibrium atmosphere having the same T e ff and logg of the dynamic 
models. The two sets of spectra are compared with Rayleigh- Jeans black bodies normalized 
to the model spectral emission for A > 8 /iin. 

It is apparent from the plot that the model spectra of ( Gem strongly depart from black 
body spectra at IRAC wavelengths: we should thus expect IRAC colors that are significantly 
different from zero magnitudes. The second most obvious characteristic of the models is the 
presence of a broad spectral feature between 4 and 6 /im, due to CO molecular absorption. 
This feature falls within the band-passes of the IRAC 4.5 and 5.8 fim filters, and is strong 
enough to have a significant effect on the IRAC colors including these bands. This feature is 
present in both the dynamic and static models, and is thus a general property of stars with 
the temperature and gravity of Cepheids, rather than being a consequence of the pulsations. 
The strength of the feature appears however to be variable with phase: it is stronger when 
the stellar atmosphere is more expanded, as at maximum radius. When the atmosphere is 
compressed (as at minimum radius, or when a shock-wave is crossing the photosphere in the 
dynamic models), the feature is reduced to a minimum. The absorption is also generally 
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stronger in the static than in the dynamic model spectra (the former corresponding to an 
atmosphere with a more "relaxed" structure, as opposed to the "compressed" atmosphere in 
the dynamic models). This is not the only difference between the static and dynamic mod- 
els: at 4>l — 0.63 (when a shock-wave is transiting in the continuum-forming photosphere) 
the dynamic and static atmospheres have very different continuum slopes in their spectral 
emission, reflecting the dramatic consequences of the pulsations in the affecting the structure 
of the Cepheid atmosphere. 

The presence of this CO absorption, and of the time-dependent hydrodynamic effects, 
induce significant color variations in the IRAC bands. Figure [7] shows the synthetic IRAC 
colors of the ( Gem dynamic and static models. The largest spread is in the [3.6] — [4.5] 
color (~ 0.8 mag and ~ 1 mag in the dynamic and static model respectively), while the 
[3.6] — [8.0] color shows a much smaller variation (~ 0.08 mag and ~ 0.02 mag). The reason 
for this difference is that the [3.6] — [4.5] color is affected by the variable CO absorption, 
while the [3.6] and [8.0] bands are outside this spectral feature. The dynamic model colors 
in fact are confined for most phases in a smaller region (~ 0.03 mag in both colors), with 
the exception of a large excursion when the shock-wave crossing the ( Gem photosphere at 
4>l = 0.63 suddenly expands the atmosphere. The [3.6] — [4.5] colors of the static models 
have instead a large variation that follows closely the changes in the effective temperature 
of the star. 

The figure also shows the actual colors of ( Gem as observed with IRAC. Due to the 
large uncertainty of this star's photometry, the measured colors cannot be used to assess the 
respective accuracy of the two class of models (the measurement is within la of both model 
colors at the observation's phase 4>l = 0.85). 

Figure [8] shows the IRAC and MIPS color-color diagrams of the program stars (including 
the 3 non- variable supergiants and 1 red giant added to the sample for comparison). The 
dashed and dotted boxes indicate the location of the dynamic and static model IRAC colors. 
It is clear that the color variation of the ( Gem models are sufficient to explain the color 
spread present for most of the program stars. The dynamic models are actually closer in 
both axes to the actual colors of the observed stars, while the static models seem to have 
a systematic offset of ~ 0.04 mag in the [3.6] — [8.0] color (the stars appear to be bluer). 
While it is unwise to extract general considerations from this effect (the color difference is 
barely significant given the photometric uncertainty, and the models are specific to only one 
of the stars), this discrepancy suggests that the color effects induced by the dynamics of 
the Cepheid pulsations may be important in determining their actual IRAC colors. Three 
stars have a [3.6] — [4.5] colors (< —0.15) significantly different from the range predicted by 
both the dynamic and static models. These stars are among the longest period Cepheids in 
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our sample. This suggest that long period Cepheids may present a deeper CO absorption, 
at least at some pulsation phases, than ( Gem. Shorter period Cepheids instead show 
small [3.6] — [4.5] color scatter, suggesting less variability in the CO feature. Polaris shows 
~ 0.1 mag redder colors than all other sources or models, but this needs to be confirmed 
with higher precision photometry for this star. 

The colors of the 4 non-variable control stars are indistinguishable from the colors of 
the Cepheids, confirming that the mid-IR CO spectral feature that is responsible for the 
color shift is common in giant stars with this spectral type even when not on the Cepheid 
instability strip. 

The comparison of the observed IRAC colors with the ( Gem models spectra indicates 
that the [3.6] — [4.5] and [3.6] — [5.8] scatter in the PC plots in Figure [5] is intrinsic. The 
[3.6] — [8.0] colors are consistent with the model colors. We therefore conclude that there is 
no measurable extrinsic infrared excess in the IRAC bands, with only the possible exception 
of Polaris. Given that [3.6] — [8.0] excess would correspond to the presence of warm dust 
(Td ~ 500 K), this means that we do not find compelling evidence, for most sources, of 
circumstellar dusty wind that could be associated with currently active dusty mass loss 



(contrary to what was stated in lMarengo et al.ll2009bl where Cepheid intrinsic color variations 



were not taken into account). This result does not exclude however the presence of winds 
with low dust content. 

The [24] — [70] color does not suffer from variability-induced scatter, but the lack of 
good photometric measurements for the 70 /xm band prevents determining reliable colors for 
most of the sources. The plot shows a trend in this color for the stars in the sample, but this 
trend is most likely a consequence of having only upper limits for the 70 //m brightness of 
most stars. Only one star (again SZ Tau) has a [24] — [70] excess above its 3cr photometric 
error, which may indicate the presence of cold dust (T^ ~ 50 K) in the proximity of the star. 



7. Discussion 



The reliability of Classical Cepheids as standard candles is of paramount importance 
for astronomy. The recent advances in infrared space astronomy have shifted the focus of 
obtaining accurate P L relations to wavelen gths longer than the visibl e, where interstellar 
extinction is reduced. iMadore et al.l (j2009bl ) and iFreedman et al.l (120091 ) have demonstrated 
how PL relations obtained at IRAC wavelengths can be effectively used to measure the 
distance of nearby galaxies. This work is a further step in the characterization of PL relations 
in all Spitzer photometric bands, including two MIPS bands at 24 and 70 /xm, by using a 
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sample of Classical Cepheids in the Galaxy. 



In order to provide a detailed comparison between th eory and ob s ervat io ns we adopted 



the la rge set o f nonlinear, connective m odels computed by lBono et all (119991 ); iMarconi et al. 



( 120051 ) and by iFiorentino et all (120071 ). For Galactic Cepheids we chose a scaled-solar chem- 
ical composition (helium, Y = 0.28; metals, Z = 0.02) and accounted for fundamental mode 
pulsators. Moreover, we covered a broad range of stellar masses (3.5 < M /M^ < H-5), and 
to account for current uncertainties affecting the size of the helium core (IBono et al.ll2006l ) 
we adopted two different mass-luminosity relations based on canonical and non-canonical 
(the latter incorporating a range of main sequence core convective overshoot) evolutionary 
models. Theoretical predictions were transfo rmed into the observation al plane using scaled- 
solar atmospheres based on ATLAS9 models (jCastelli &: Kuruczll2003l ). and multiplied with 
the IRAC band-passes. The slopes and zero points of these model Leavitt laws are listed 
in the last two columns of Table [61 While the zero points are in general agreement with all 
our fits, the slopes are significantly shallower than our best fits obtained with both "new" 
and "old" distances (the disagreement is however reduced in the fits adopting the "new" dis- 
tances). The PL relations derived with the exclusive use of the astrometric distances (that 
do not rely on the p(P) relation), are however in excellent agreement with the theoretical 
relations. This result may suggest that the period dependence of the p— factor currently 
adopted in IRSB distances may still need further refinement. Because of the small number 
of Cepheids and the limited range of periods for which accurate parallaxes are available, we 
could not attempt to invert the problem and estimate the p(P) dependence with our data. 

Comparison between our Galactic PL relation s and the relations derived for the LMC 
by iMadore et all (j2009al ) and iNgeow et all (120091 ) leads to contradictory results. The PL 
relations obtained with both the "new" and "old" IRSB distances are steeper than the 
LMC relations, while the Galactic PL relations we obtain with the astrometric distances are 
sig nificantly sha ll ower. To a dd to the conf u sion, it should be noted that the slopes derived 
by INgeow et all (120091 ) and IMadore et all (j2009al ) disagree by more than their respective 
uncertainties, despite using stars from the same galax}o The contradictory results in our 
Galactic PL relation slopes shows once more the effects of the systematics introduced by 
the IRSB distances and their dependence on the p(P) relation. Based on the more reliable 
astrometric distances we should conclude that the slope of the PL relation could be shallower 
at Galactic metallicity than in the LMC, but even this result is not statistically significant, 
due to the larger uncertainty in the PL slope. In conclusion, the insufficient reliability of the 



3 T he difference between the two LMC results can be an issue of crowding and binarity in the lNgeow et al. 
(|2009| ) sample (selected without individual image inspection, and containing a larger fraction of fainter short 
period Cepheids) that could result in shallower slopes. 
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IRSB distances, and the small number of stars for which astrometric distances are available, 
prevents us from resolving the dependence of the PL relations from metallicity. This problem 
has already been noticed at optical and near-IR wavelengt h s by other autho rs (see e.g. 
Storm et alihooi iFouque et alJl2007t iRomaniello et al.ll200a l. iRiess et all J2009h also noted 
that the precision of current datasets does not allow measurement of the dependence of 
the PL relation on metallicity; those authors measured statistically consistent slopes at 1.6 
micron for the Milky Way, the LMC and other galaxies. 

We also cannot determine a reliable wavelength dependence of the PL slope, as all our 
values are within la from the average slope value. This is consistent with the expected 
flattening of the PL slope wavelength depende nce in the mid-IR sho wn by the models. If 
we compare our slope values with Figure 4 in iFreedman et al.1 (120081 ). however, we note a 
similar trend. Our result for 4.5 /im is marginally higher (by ~ 0.1) than the slope at the 3.6 
and 8.0 /im, as also shown by the LMC Cepheid fit (and by the models). This anomaly in 
the 4.5 /im (and, to a minor extent, 5.8 /im) slope is related to the presence of the variable 
CO band, which increases the amplitude variation in these bandpasses. 

For the first time we have measured the PL relation in the MIPS 24 and 70 /im bands. 
The 24 /im relation, in particular, is important for the determination of C epheid distances 
in the Galactic plane from datasets such as MIPSGAL (jCarey et al.ll2009l ). where high ex- 
tinction from ISM dust may complicate their detection at IRAC bands. 

The determination of a PC relation is in principle not limited by the uncertainty in the 
distances affecting our PL relations. Figure [5j however, clearly shows that such a relation 
can be derived with sufficient accuracy only in the [3.6] — [8.0] color. Other colors show 
a scatter that is increasingly larger with the period, as high as ~ 0.1 mag, one order of 
magnitud e larger than our c olor accuracy of ~ 0.02. The fact that this scatter was not 
found by iNgeow et al.1 (120091 ) may again be a consequence of the difference in the median 
period between our Galactic and their LMC sample. 

Based on our data and on numerical modeling of one representative Cepheid (£ Gem), 
we conclude that the color scatter affecting long period Cepheids is a consequence of variable 
stellar CO absorption, rather than infrared excess from circumstellar dust. Circumstellar 
dust at ~ 500 K (the temperature where dust thermal emission is maximum in the IRAC 
bands) would produce the largest excess in the [3.6] — [8.0] color, which we do not observe, 
with the possible exception of a weak detection (~ 0.1 mag) around Polaris. Our ( Gem 
models show that scatter in the IRAC 4.5 and 5.8 //.m bands can be induced by variations in 
effective temperature and/or the propagation of shocks through the Cepheid atmospheres. A 
more comprehensive modeling effort, including Cepheids with different periods and pulsation 
modes, is however required to confirm and quantify this effect. 
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Not finding an infrared excess in most, if not all, stars is somewhat une xpected, given 



the de tection of circumstellar shells by near-IR interferometers. Figure 8 in IMerand et al. 



(120071 ) shows a number of targets, including the long period Cepheids Y Oph and t Car 
(both part of our sample), with a significant (4-5%) excess in the K band. We do not 
measure any excess of this magnitude in these or other stars, which seems to rule out the 
presence of circumstellar dust in the shells detected around these and other Cepheids. The 
interferometric determination of t he radius of these ci rcumstellar shells, however, indicate 
values as low as ~ 2 stellar radii (IMerand et al.ll2006l ). For a Cepheid with T g g 
6,000 K, dust equilibrium temperature at ~ 2 R* would be Td — T* 



5,000- 



2 R* would be T d ~ T*/^R d /R* < 3,500 K. 
This is well above the sublimation temperature of any known astronomical dust. 

Our results can be reconciled with the interferometric detection of circumstellar shells 
close to the star if the emission is not due to dust, but rather t o som e st rong molecular 



l ine em ission. Given that the observations in IMerand et all ( 120061 . 120071 ) and iKervella et al. 



( 120061 ) were made in the K band, a possible candidate for this emission is shocked H 2 . The 
presence of H 2 shocked emission lines could further co ntribute to the sca tter in the IRAC 
4.5 fim band, for certain densities and shock velocities (jSmith et al.ll2006l ). This hypothesis 
supports the idea that Cepheid stars indeed have a str ong stellar win d asso ciated with a 
pulsation driven mass-loss mechanism, as suggested by IMerand et al.1 (120071 ). Unlike the 
case of mass loss in red giants, AGB and supergiant stars (showing a dust-driven wind), 
this stellar wind would be largely dustless, due to the higher temperature of the star. This 
hypothesis needs to be tested by near-IR multi-epoch spectral monitoring of these stars. 

This hypothesis does not imply that a Cepheid stellar wind is completely devoid of dust, 
that may be condensing at larger radii, in quantities below the detection limits of our IRAC 
observations (see I Car mid-IR detection of mid-IR circumstellar emission with the VLT 



MIDI and VISIR instruments, IKervella et al.ll2009i ). Cold dust may also be collected by the 
outflow, at much larger distances ( thousands of AU), fr om the ISM, as in the case of the 
well known nebula around RS Pup (IKervella et al.ll2009i ). We find evidence supporting this 
mechanism in our IRAC and MIPS images. Extended emission at 5.8, 8.0, 24 and 70 \xm 
around S Cep indicate that this star may be losing mass due to a strong wind pushing into the 
local interstellar medium, which is leading to the formation of a 70 /mi bow shock detected 
at large distance from the star (~ 10,000 AU). A detailed analysis of this phenomenon is 
being published elsew here (Marengo et al., in preparation). In a more comprehensive paper 
(fearmbv et al]|200gh we will discuss the presence of spatially resolved extended emission at 
24 and 70 /xm around more targets, and quantify the occurrence of mass loss in the Cepheid 
phase that can be inferred from our imaging data. 
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8. Summary 

We have derived the PL and PC relations in Spitzer /IRAC bands for a sample of 
Galactic Cepheids. These relations are critically dependent on the choice of the period 
dependence of the p-factors used for the distance determination with the IRSB method, 
even though the uncertainties in our fits prevent the assessment of the best p(P) relation 
choice. The best agreement with theoretical PL relations, however, is obtained when only 
distances obtained by astrometric methods are used. We do not detect statistically significant 
variations between the slope and zero points of t he PL relations between our Galactic sample 



and LMC relations obtained by iMadore et al.l (j2009al ). despite the difference in metallicity. 



We find that the intrinsic variations in the 4.5 and 5.8 jum fluxes are larger for long period 
Cepheids. These variations (of the order of ~ 0.1 mag) are related to deep CO absorption, 
dependent on the stellar T e ff and hydrodynamic effects associated to the stellar pulsations. 
We do not find significant infrared excess related to warm circumstellar dust, except for a 
weak excess detected at IRAC wavelengths for Polaris, and at 70 fim for SZ Tau. This may 
rule out the presence of extensive dust driven mass loss in the Cepheid phase, but leaves 
open the possibility of pulsation-driven mass loss from a dust-poor wind, as suggested by 
recent interferometric observations. 
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This preprint was prepared with the AAS IATgX macros v5.2. 



Table 1. Target list and catalog of observations 



Name 


IRAC sub JD a 


IRAC full JD a 


MIPS JD a 




D [kpc] b 




D [kpc] c 


P[d] 


Mode 


DT Cyg 


54065 


535 




53935.208 





457±0 


069 





535±0 


029 d 


2.54 


1st 


SZ Tau 


54006 


076 


54189.135 


54400.742 





513±0 


029 





556±0 


034 c 


3.15 


1st 


RT Aur 


54036 


201 




54202.804 





417±0 


033 f 





417±0 


033 f 


3.73 


F 


Polaris 1 






52991.065 


54024.680 





130±0 


002 f 





130±0 


002 f 


3.97 


1st 


BF Oph 


53997 


204 




53979.664 





637±0 


032 





714±0 


011 


4.07 


F 


FF AqP 


54039 


420 




54023.113 





356±0 


023 f 





356±0 


023 f 


4.47 


F 


V350 Sgr 1 


54002 


135 




54013.655 





935±0 


052 





979±0 


054 d 


5.15 


F 


5 Cep 






53957.101 


53935.018 





273±0 


01 l f 





273±0 


01 l f 


5.37 


F 


V Cen 


53957 


367 


53957.365 


53980.875 





606±0 


044 





684±0 


020 


5.49 


F 


BB Sgr 


54002 


629 




54013.643 





787±0 


031 





801±0 


010 


6.64 


F 


U Sgr 


53999 


677 


53997.790 


54009.048 





565±0 


019 





585±0 


006 


6.75 


F 


V636 Sco 1 


53997 


194 




53983.521 











839±0 


046 d 


6.80 


F 


U AqP 


54039 


442 




54016.545 





606±0 


029 





691±0 


038 d 


7.02 


F 


f] Aql 






54039.439 


54238.612 





241±0 


014 





250±0 


005 


7.18 


F 


W Sgr 1 


53997 


264 




54013.542 





439±0 


038 f 





439±0 


038 f 


7.59 


F 


GH Lup 


53961 


333 




53979.733 


1 


124±0 


203 


1 


023±0 


051 s 


9.28 


F 


S Mus 1 


53960 


571 




53935.332 





820±0 


054 





916±0 


046 s 


9.66 


F 


S Nor 


53999 


638 


54351.249 


53983.802 





943±0 


044 





959±0 


014 


9.75 


F 


Dor 






53960.988 


53935.322 





318±0 


016 f 





318±0 


016 f 


9.84 


F 


C Gem 






54228.830 


54046.15 





360±0 


023 f 





360±0 


023 f 


10.15 


F 


XCyg 


54065 


539 




53935.023 


1 


163±0 


027 


1 


213±0 


010 


16.39 


F 


Y Oph 


53997 


328 




54202.088 





552±0 


040 





573±0 


008 h 


17.13 


F 


VY Car 


53957 


067 




53935.360 


1 


818±0 


099 


1 


994±0 


020 


18.91 


F 


SW Vol 


54098 


538 




54202.474 


2 


381±0 


057 


2 


506±0 


029 


23.44 


F 


T Mon 1 


54395 


313 




54021.550 


1 


389±0 


058 


1 


455±0 


037 


27.02 


F 


AQ Pup 


54228 


851 




54107.409 


3 


030±0 


184 


3 


194±0 


066 


30.10 


F 


£ Car 






53960.984 


53935.328 





498±0 


050 f 





498±0 


050 f 


35.55 


F 


U Car 


53960 


567 




53935.349 


1 


492±0 


067 


1 


565±0 


023 


38.77 


F 



Table 1 — Continued 



Name 


IRAC sub JD a 


IRAC full JD a 


MIPS JD a 


D [kpc] b 


D [kpc] c 


P[d] 


Mode 


RS Pup 


54098.477 




54073.425 


1.818±0.099 


2.052±0.151 


41.39 


F 


HD183864J 


54064.419 




54046.098 










HD182296 j 


54226.294 




54018.537 










7 Phe* 


53957.014 




53935.315 










ip And j 




54094.302 


54106.889 











a JD-2400000.5 at observation start 



b Distance from lFouque et al.l (|2007l ) unless noted 



c Distance from lStorm et al.l (|2004h unless noted 
d This work (Table© 



°Distance from lGroenewegen et al.l (|2004r ) 



f Distance from iBenedict et al.l (|2007l ) or Ivan Leeuwen et al.l (|2007l ) for Polaris 

s Distance from iRomaniello et al.l (|2008l ) (with 5% assumed uncertainty) 

h Distance from Barnes et al. ( 20051) . Least Square fit 

'Spectroscopic binary 

^Supergiant or M giant (7 Phe) star 



Table 2. "Old" distances missing from previous catalogs 



Name 



Ro[Rq 



8q [mas] 



V 



K 



E(B - V) 



V 



K 



d [kpc] 



DT Cyg 30.7±0.9 

V350 Sgr 40.7±4.5 

V636 Sco 50.0±4.8 

U Aql 51.3±6.0 



0.570±0.018 
0.424±0.047 
0.568±0.054 
0.830±0.097 



5.744±0.010 
7.483±0.010 
6.654±0.010 
6.446±0.010 



4.430±0.010 
5.141±0.010 
4.409±0.010 
3.893±0.010 



0.039±0.020 
0.312±0.020 
0.217±0.020 
0.399±0.020 



5.635±0.010 
6.516±0.010 
5.981±0.010 
5.209±0.010 



4.416±0.010 
5.031±0.010 
4.332±0.010 
3.752±0.010 



0.535±0.029 
0.979±0.054 
0.839±0.046 
0.691±0.038 



Table 3. IRAC and MIPS photometric calibration 



Item 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


Isophotal A [pm] 


3.550 


4.493 


5.731 


7.872 


23.68 


71.42 


FLUXCONV [(MJy/sr)/(DN/s)] 


0.1088 


0.1388 


0.5952 


0.2021 


0.0447 


702. 


6.1" aperture corr. [mag] a 


0.051±0.002 


0.055±0.002 


0.056±0.006 


0.065±0.005 






fcorr (subarray) 


0.994 


0.995 


0.939 


1.006 






F„(Vega) [Jy] 


280.9±4.1 


179.7±2.6 


115.0±1.7 


64.1±0.9 


7.14±0.08 


0.775±0.009 



a Aperture used for FF Aql, VY Car and AQ Pup IRAC subarray data 
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Table 4. IRAC aperture vs. PSF fitting photometry comparison 3. 



Name 


IRAC mode 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


V Cen 


subarray 


4.336±0.010 


4.356±0.010 


4.377±0.003 


4.347±0.003 




full frame 


4.332±0.270 


4.392±0.156 


4.392±0.093 


4.362±0.060 


U Sgr 


subarray 


3.811±0.010 


3.856±0.010 


3.864±0.008 


3.829±0.006 




full frame 


3.844±0.236 


3.772±0.105 


3.789±0.054 


3.789±0.053 


S Nor 


subarray 


3.969±0.010 


4.015±0.010 


4.035±0.008 


4.000±0.006 




full frame 


3.984±0.150 


4.072±0.136 


4.027±0.067 


4.005±0.044 


SZ Tau 


subarray 


4.165±0.010 


4.164±0.010 


4.206±0.009 


4.185±0.007 




full frame 


4.144±0.200 


4.144±0.200 


4.221±0.080 


4.195±0.076 


Avg diff. 




-0.007±0.021 


0.009±0.065 


0.021±0.049 


0.003±0.022 



a With the exception of V Ccn, full frame and subarray data of each source have been 
obtained at different epochs, and the difference in photometry can be related not only on 
photometry biases, but also to the source variability 



Table 5. Final Spitzer photometry of all targets 



Name 


mode 




\"i fil 






[4.5J 






\k el 
[5.8J 






[8.UJ 






[24j 




[7UJ 




DT Cyg 


subarray 


4. 


,340±0 


.010 


4 


.342±0. 


.011 


4 


.381±0. 


.010 


4. 


,386±0 


.008 


4, 


,367±0 


.002 


3.791±0 


.257 


SZ Tau 


subarray 


4. 


165±0 


.010 


4 


,164±0. 


.010 


4 


.206±0. 


.009 


4. 


185±0 


.007 


4. 


209±0 


.002 


3.013±0 


.195 


RT Aur 


subarray 


3. 


.865±0 


.010 


3 


.880±0. 


.010 


3 


,924±0. 


.008 


3. 


,914±0 


.006 


3. 


773±0 


.002 


3 


,645 


Polaris 


full frame 





,573±0 


.037 





.493±0 


.051 





,411±0. 


.070 


0. 


,459±0 


.042 


0. 


,723±0 


.001 


0.506±0 


.050 


BF Oph 


subarray 


5. 


141±0 


.010 


5 


,136±0. 


.010 


5 


,171±0. 


.005 


5. 


162±0 


.004 


5, 


,031±0 


.005 


3 


.698 


FF Aql 


subarray 


3. 


.391±0 


.010 


3 


.390±0. 


.010 


3 


.436±0. 


.008 


3. 


,428±0 


.007 


3, 


,338±0 


.001 


3, 


.458 


V350 Sgr 


subarrav 


4 


,948±0 


.010 


4 


.957±0 


.010 


4 


.988±0. 


.004 


4. 


,988±0 


.003 


5, 


,032±0 


.009 


3, 


.991 


5 Cep 


full frame 


2. 


174±0 


.044 


2 


,182±0. 


.037 


2 


.166±0. 


.032 


2. 


150±0 


.039 


2. 


120±0 


.001 


2. 


.174 


V Cen 


subarray 


4, 


,336±0 


.010 


4. 


.356±0 


.010 


4 


,377±0. 


.003 


4. 


,347±0 


.003 


4, 


,424±0 


.003 


2. 


.813 


BB Sgr 


subarrav 


4. 


350±0 


.010 


4 


.380±0. 


.010 


4 


.386±0 


.003 


4. 


370±0 


.003 


4. 


339±0 


.004 


3. 


.814 


U Sgr 


subarray 


3 


,811±0 


.010 


3 


.856±0 


.010 


3 


.864±0. 


.008 


3 


829±0 


.006 


3. 


760±0 


.003 


1. 


.910 


V636 Sco 


subarray 


4. 


,375±0 


.010 


4 


.389±0. 


.011 


4 


.434±0. 


.010 


4. 


400±0 


.008 


4. 


384±0 


.001 


2. 


.913 


U Aql 


subarray 


3. 


.810±0 


.010 


3 


,869±0 


.010 


3 


,877±0. 


.008 


3. 


839±0 


.006 


3. 


665±0 


.002 


3.682±0 


.294 


n Aql 


full frame 


2. 


,003±0 


.041 


1 


,997±0. 


.048 


1 


990±0 


.021 


1 


.970±0 


.047 


1. 


856±0 


.001 


1.860±0 


.165 


W Sgr 


subarray 


2. 


.787±0 


.010 


2 


,855±0. 


.010 


2 


.877±0. 


.005 


2. 


825±0 


.004 


2. 


808±0 


.001 


2. 


.154 


GH Lup 


subarray 


4. 


,692±0 


.010 


4. 


.716±0. 


.010 


4 


,731±0. 


.004 


4. 


706±0 


.003 


4. 


,644±0 


.006 


2. 


.085 


S Mus 


subarray 


3 


,835±0 


.010 


3 


,845±0 


.010 


3 


,884±0. 


.008 


3 


.877±0 


.006 


3. 


902±0 


.001 


3.156±0 


.278 


S Nor 


o LI 1 Ki'l i. (X y 


3. 


,969±0 


.010 


4 


015±0. 


.010 


4 


035±0. 


.008 


4 


ooo±o 


.006 


4 


086±0 


.002 


3 


.102 


13 Dor 


full frame 


1 


,899±0 


.059 


1 


,990±0 


.068 


1 


,937±0. 


.049 


1 


.880±0 


.037 


1. 


.858±0 


.001 


1.864±0 


.103 


C Gem 


full frame 


2. 


,045±0 


.039 


2 


,045±0. 


.032 


2 


.060±0. 


.033 


2. 


045±0 


.032 


1. 


982±0 


.001 


1.799±0 


.124 


XCyg 


subarray 


3 


,741±0 


.010 


3 


,880±0. 


.010 


3 


,841±0. 


.007 


3 


.765±0 


.006 


3. 


680±0 


.001 


3 


.352 


Y Oph 


subarray 


2. 


,543±0 


.010 


2 


,510±0. 


.010 


2 


.566±0. 


.004 


2. 


556±0 


.003 


2. 


509±0 


.001 


2. 


.965 


VY Car 


subarray 


4. 


,497±0 


.010 


4 


,591±0. 


.010 


4 


.580±0. 


.007 


4. 


505±0 


.006 


4. 


527±0 


.006 


2. 


.391 


SW Vcl 


subarray 


4. 


,891±0 


.010 


4 


,908±0. 


.010 


4 


,919±0. 


.006 


4. 


.887±0 


.003 


5. 


,064±0 


.003 


4. 


.302 


T Mon 


subarray 


3 


,309±0 


.010 


3 


.469±0. 


.010 


3 


,428±0. 


.006 


3. 


320±0 


.005 


3. 


,422±0 


.001 


2. 


.978 


AQ Pup 


subarray 


5. 


,014±0 


.010 


4 


.982±0 


.010 


5. 


.009±0. 


.007 


4. 


.993±0 


.006 


4. 


.889±0 


.004 


4. 


.171 


£ Car 


full frame 


1. 


.020±0 


.028 


1 


020±0 


.049 


1 


.003±0 


.048 


0. 


927±0 


.033 


0. 


,720±0 


.001 


0.690±0 


.053 



Table 5 — Continued 



Name 


IRAC mode 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


U Car 


subarray 


3.543±0.010 


3.694±0.010 


3.654±0.007 


3.535±0.005 


3.187±0.006 


0.454 


RS Pup 


subarray 


3.324±0.010 


3.336±0.010 


3.350±0.006 


3.331±0.005 


3.316±0.001 


-0.068 


HD183864 


subarray 


4.281±0.010 


4.311±0.011 


4.321±0.010 


4.310±0.007 


4.247±0.002 


3.959 


HD182296 


subarray 


4.006±0.010 


4.076±0.010 


4.098±0.009 


4.038±0.006 


3.993±0.002 


3.361 


7 Phe 


subarray 


-0.627±0.048 


-0.526±0.149 


-0.504±0.121 


-0.531±0.211 






ip And 


full frame 


2.415±0.076 


2.507±0.084 


2.462±0.037 


2.389±0.034 


2.332±0.001 


2.307±0.157 



Note. - 



Magnitudes without errors are lower limits (flux upper limits) 



Table 6. PL relation coefficients for different samples 



■'New" distances 3. 



Astr. distances' 3 



''Old" distances 



Model (Z=0.02, Y=0.28) c 



Band 


a 


P 


Ct 


a 


a. 


13 


a 


13 


3.6 


-3.47 ±0.06 


-5.72 ±0.07 


— 3.16 ± 0.22 


-5.74 ±0.18 


— 3.54 ± 0.04 


-5.76 ±0.04 


-3.13 ±0.01 


-5.79 ±0.01 


4.5 


-3.38 ±0.06 


-5.68 ±0.07 


— 3.06 ± 0.23 


-5.74 ±0.19 


-3.45 ±0.04 


-5.72 ±0.04 


-3.04 ±0.01 


-5.71 ±0.01 


5.8 


-3.44 ±0.06 


-5.68 ±0.07 


-3.10 ±0.23 


-5.75 ±0.19 


-3.51 ±0.03 


-5.70 ±0.04 


-3.07 ±0.01 


-5.76 ±0.01 


8.0 


-3.46 ±0.06 


-5.73 ±0.07 


-3.16 ±0.22 


-5.80 ±0.18 


-3.57 ±0.03 


-5.74 ±0.04 


-3.12 ±0.01 


-5.80 ±0.01 


24 


-3.52 ±0.06 


-5.71 ±0.06 


-3.51 ±0.21 


-5.75 ±0.17 


-3.67 ±0.03 


-5.78 ±0.03 






70 


-3.18 ±0.26 


-5.87 ±0.22 


-3.34 ±0.27 


-5.87 ±0.23 


-3.10 ±0.25 


-5.89 ±0.22 






K 


-3.37 ±0.06 


-5.65 ±0.02 


-3.32 ±0.12 


-5.71 ±0.03 


-3.67 ±0.12 


-5.69 ±0.03 







a Distances from lFouque et al. ( 2007 ). Benedict et al. ( 2007 ) or van Leeuwen et al. ( 2007 ) 



Distances from iBenedict et al.1 ( 20071 ) or Ivan Leeuwen et al.l ( 20071 ) 



Distances from IStorm et al.l (|2004[ ). iGroenewegen et alj (|2004l ). iBarnes et al.l (|2005l ), IBenedict et al.l (|2007l) . Ivan Leeuwen et al 



(|2007f ) or this work (Table E]) 

d Based on the models described in Section [7] 



DT Cyg SZ Tau RT Aur BF Oph FF Aql 




AQ Pup U Car RS Pup 



WWW 



Fig. 1. — Thumbnail images of the 23 Cepheids observed in the IRAC subarray mode 
at 3.6 jum, shown with a logarithmic color scale. The size of each thumbnail is about 
45 x 45 arcsec. Only three targets (FF Aql, VY Car and AQ Pup) show one or more 
background stars within the standard 12.2 arcsec radius aperture used to measure the stars 
photometry. For these three stars a smaller 6.1 arcsec radius aperture has been used. 
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Fig. 2. — Comparison between aperture and PSF fitting photometry of four stars in our 
sample: #1 V Cen, #2 U Sgr, #3 S Nor and #4 SZ Tau. The error bars are the combined 
uncertainty of the offset between the aperture and PSF-fitting magnitudes. The dashed 
line is the average offset between the two sets of measurements (much smaller than the 
photometric error of each data-point). Note that only V Cen was observed in both subarray 
and full frame at the same epoch. 
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0.5 




log P [days] 



Fig. 3. — Period vs. absolute magnitude plots for all IRAC and MIPS bands. Circles are 
fundamental mode Cepheids while triangles are first o vertone (with period "fundamental- 
ized" using the relation from iFeast fc Catchpold 119971 ). This plot shows 3a error bars for 
IRAC subarray photometry and the uncertainty listed in Table [T] for all other datapoints . 



The solid line shows the best-fit PL relation using all the distances from lFouque et al 
while the dashed line uses only stars with astrometric distances. 



2003), 
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AM[8.0] AM[3. 6-8.0] 

Fig. 4. — Left panel — Correlated residuals from the PL relation. Right panel — Color 
residuals from the PL relation. Circles are fundamental mode Cepheids, while triangles are 
first overtone. The two dashed lines show the 3<ta scatter around the unity-slope line (solid 
line) . 
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Fig. 5. — Spitzer colors vs. period for our sample. The solid line in the [3.6] — [4.5] color 
panel is the best fit PC relation. The solid line in the other panels, when present, is a marker 
for the zero color of Rayleigh- Jeans atmospheres). The dashed and dotted lines indicate the 
1 and 3 a dispersion expected for the 3.6 /im lightcurve amplitude estimated in Section loTTl 
Circles are fundamental mode Cepheids and triangles first overtone. SZ Tau has a [24] — [70] 
color excess above 3<r the noise level, and Polaris has consistently red (~ 0.1 mag) IRAC 
colors. 
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Fig. 6. — Mid-IR model spectra of the 10 d perio d Cepheid C Gem. Th e thick solid lines are 
time-dependent hydrodynamic models derived by lMarengo et al.l (120021 ). The thin solid lines 
are instead equivalent models having the same T e ff and logg as the dynamic models, but 
an atmosphere in hydrostatic equilibrium. The dotted lines are Rayleigh- Jeans blackbody 
fits of the spectra for A > 8 //m. The spectra are plotted for a number of significant phases: 
maximum V luminosity (0x = 0), maximum radius (0^ = 0.32), minimum V luminosity 
(0l — 0.49), phase at which a shock- wave is crossing the photosphere (0^ = 0.63), minimum 
radius (0^ = 0.74) and the phase at which the star was observed with IRAC (0^ = 0.85). A 
broad CO spectral feature is present at most phases, but with different strength, within the 
4.5 and 5.8 /zm IRAC pass-bands. 
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Fig. 7. — IRAC colors of the dynamic (solid line) and static (dotted line) models of £ Gem. 
The jitter on the model tracks is due to the finite resolution of the T e ff and log g model grid. 
The colors of the dynamic model at the same phases of the spectra plotted in Figure [6] are 
marked (square points). The model and actual colors of the source at the epoch of the IRAC 
observations are also plotted. The dashed and dotted square boxes indicate the "bounding 
boxes" of the two class of models (dynamic and static respectively) in the IRAC color space. 
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-0.1 -0.05 0.05 0.1 -0.1 -0.05 0.05 0.1 



[3.6]-[8.0] [3.6]-[8.0] 

Fig. 8. — Color-color diagrams of our sample. Left panel — IRAC [3.6] — [8.0] vs. [3.4] — [4.5] 
colors. The [3.6] — [8.0] color shows a small spread (~ 0.05 mag) for most of the sources. 
The [3.6] — [4.5] color has instead a much larger spread, well above the variations allowed 
by the photometric uncertainties. Polaris is the source with redder colors (~ 0.1 mag). The 
dashed and dotted boxes indicate the location of the ( Gem dynamic and static models, 
respectively. Right panel — IRAC [3.6] — [8.0] vs. MIPS [24] — [70] color: only one source 
(SZ Tau, triangle indicated by arrow) has a statistically significant excess in the [24] — [70] 
color. Other sources seem to be on a sequence of increasing colors, but the trend may be 
a spurious artifact, due to the absence of reliable 70 /iin detections. The 3 supergiants and 
the red giant in our sample (square symbols) have indistinguishable colors from those of the 
Cepheid stars, in both diagrams. 



